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^ (57) Abstract: Disclosed are a method and a device for the differential spectrophotometry of non-clear media by means of spectro- 
^ temporal imaging in counting mode. In order to analyze a non-clear medium, said medium is illuminated by at least one light pulse, a 
spectral and temporal transmission image is obtained of the illuminated medium in counting mode, and the image and the derivatives 
O thereof are processed in order to obtain information on the non-clear medium. The invention applies to the analysis of diffusing and 
^ absorbing media, e.g. milk. 

^ (57) Abrege : Procddf et dispositif de spectrophotom6trie differennelle de milieux non limpides par imagerie spectro-temporelle en 
mode comptage. Selon l'invenuon, pour analyser un milieu non limpide, on eclaire ce milieu par au moins une impulsion lumineuse, 
Q on acquiert, a partir du milieu ainsi eclaire\ une image spectrale et temporelle de transmission, en mode comptage, et Ton traite 
£^ l'image et ses derivSes pour acqudrir des informations sur le milieu non limpide. L'invention s'applique a 1' analyse de milieux 
^ diffusants et absorbants, par exemple le lait. 



WO 2004/013617 



PCT/FR2003/002411 



PROCEDE ET DISPOSITIF DE SPECTROPHOTOMETRIE 
DIFFERENTIELLE DE MILIEUX NON LIMPIDES PAR IMAGERIE 
S PECTRO - TEMPORELLE EN MODE COMPTAGE 



DESCRIPTION 



DOMAINE TECHNIQUE 



La presente invention concerne un proc§d£ 
et un dispositif de spectrophotometrie impulsionnelle . 

L 1 invention s 1 applique & 1 1 analyse de 
10 milieux qui ne sont pas limpides et trouve des 
applications dans divers domaines, notamment : 

- I 1 analyse du lait et des produits d£rivSs 
de ce dernier , 

le contrdle de produits tels que les 
15 viandes, les ceufs, les fruits et les poissons, 

- le contr61e des boues, 

- le contrdle des dechets industriels, 

- le domaine medical et 

- 1* analyse des milieux granulaires, poreux 
20 ou fractures. 



ETAT DE LA TECHNIQUE ANTERIEURE 

On connait d§ja diverses techniques de 
spectrophotometrie de milieux limpides. 
25 On connait en particulier des techniques de 

spectrophotometrie de milieux limpides, permettant de 
mesurer l 1 absorption sur 1' ensemble d f un spectre, sans 
balayage, en particulier & l'aide de dgtecteurs 
comportant des barrettes de photodiodes. 
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Ces techniques autorisent la 

spectrophotomStrie dif f erentielle et les analyses de 
melanges de colorants en temps reel, dans le cas de 
tels milieux limpides et meme dans le cas de milieux 
5 trSs peu diffusifs. 

Precisons d£s maintenant que la prSsente 
invention permet de determiner les derivees partielles 
premidres, secondes et n-ieme des flux photoniques 
mesurSs par rapport & la longueur d'onde et au temps, 

10 . en plus des informations que ces techniques connues 
permettent d'obtenir* 

On connait aussi diverses techniques de 
spectroscopie rSsolue en temps ("time -resolved 
spectroscopy") : pour la fluorescence resolue en temps, 

15 on a dejl propose l 1 imagerie spectro-temporelle de 
fluorescence. 

Une telle, imagerie donne accds & des 
distributions spectrales et temporelles, qui sont 
specif iques de l 1 analyse des temps de d^clin ("decay 

20 times") de fluorescence & plusieurs composantes. 

II convient cependant de noter que 
1'intSrit de la spectrophotom§trie impulsionnelle 
dif f erentielle pour 1 1 imagerie spectro-temporelle de 
transmittance, avec une source lumineuse impulsionnelle 

25 a large spectre, n f est ni divulgug ni suggSre dans 
l'art anterieur (tel que le document WO 01/22063 Al) • 
Seules les resolutions spectrales et temporelles 
avaient et6 explicitees mais jamais 1' utilisation des 
dSrivees partielles en fonction du temps de vol et/ou 

30 du spectre n'avait 6t€ envisagee en vue d'une analyse 
chimique et physique des milieux non limpides. 
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II convient en outre de noter que, dans le 



cadre de l ! imagerie spectro-temporelle de transmittance 
des milieux non limpides, le mode de comptage de photo- 
electrons uniques ("single photo-electrons") , ou SPE, 
5 est aussi puissant que pour l'imagerie spectro- 
temporelle de fluorescence car il permet de mesurer les 
flux photoniques non en uniti arbitraire mais en unit6 
de comptage d'evSnement. Ceci permet aussi une plus 
grande dynamique de mesure. 

10 Le mode de comptage est r§alise par 

binarisation de 1' image par seuillage. Cette 
binarisation provoque la mise & zero du bruit de la 
chalne de detection. Cette operation non lin^aire de 
seuillage autorise 1' utilisation d'une fenetre 

15 d' integration sur tin temps trds long (pouvant dgpasser 
plusieurs heures) . Ce mode permet de compter des 
evenements rares lies aux photons non absorbes ayant 
sSjournS un long moment dans le milieu d'gtude. Ces 
photons « rescapSs » portent une tr£s forte valeur 

20 informative sur le volume sonde. La detection d'une 
zone de pixel attribuable £ ce photo-electron est 
realisee par seuillage et traitement adequat. Le mode 
comptage incorpore la reduction de cette zone (en forme 
de disque ou de geometric specifique) £ un seul pixel 

25 ou k des coordonn€es a une Schelle subpixel. 

On se reportera aussi au document suivant : 
S. Anderson-Engels et al., "Multispectral 
tissue characterization with time resolved detection of 
30 diffusely scattered white light" Optics Letters (1993) 
18, pp. 1697-1699. 
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Cependant, ce document a £t£ trds peu 
exploite et reste limite au domaine de I'optique des 
tissus ("tissue optics"), sans aucune mention d'une 
application a la chimie analytique ni d' aspects 
5 probabilistes et statistiques de ce type d'imagerie 
(tels que le comptage de SPE et 1' utilisation des 
d6riv€es partielles) . 

Pour ce qui concerne le comptage, on se 
reportera aussi au document suivant et aux brevets 
10 d6pos£s par Hamamatsu au dSbut des ann£es 90: 

Watanabe M. , Koishi M. , Fujiwara M. , 
Takeshita T. and Cieslik W. (1994) Development of a new 
fluorescence decay measurement system using two- 
dimensional single-photon counting. J". Photochem. 
15 Photobiol. AzChem 80,. 429-432. 

Cependant ce document est orient^ vers le 
« two-dimensional single-photon counting » en emission 
de lumiere par fluorescence ou en tout autre emission a 
une autre longueur d'onde que la source initiale mais 
20 jamais en transmittance et en particulier en 
transmittance d'un laser blanc impulsionnel. 



EXPOSE DE L' INVENTION 

La presente invention propose un proced& et 
25 un dispositif d' analyse d'un milieu non limpide. 

En particulier, 1 ' invention vise k mesurer 
globalement et in situ des concentrations, ou des 
variations de concentrations, d'absorbeurs et de 
diffuseurs dans un milieu non limpide. 
30 De fa<?on precise, la present e invention a 

pour objet un proc£d£ de spectrophotometrie 
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dif fSrenti ell e, pour analyser un milieu non limpide, ce 
proc£d§ etant caracterise en ce que : 

- on iclaire le milieu non limpide par au 
moins une impulsion lumineuse permettant 1 ' utilisation 
5 ultSrieure de la d6rivee partielle, en fonction de la 
longueur d'onde, d'au moins une image spectro- 
temporelle de transmission acquise & partir du milieu 
ainsi 6claire, 

on acquiert, Sl travers au moins tin 

10 collecteur de lumi^re, k partir du milieu ainsi 
6clair§, au moins une image spectro-temporelle de 
transmission en mode comptage permettant 1' utilisation 
ultSrieure des derivees partielles de 1' image en 
fonction de la longueur d'onde et du temps de vol de 

15 1' impulsion lumineuse, et 

on traite 1' image et ses derivees 
partielles en fonction de la longueur d'onde et du 
temps de vol pour acqu§rir des informations sur le 
milieu non limpide. 

20 Selon un mode de mise en oeuvre prefSr§ du 

proc£dS objet de 1' invention, on exploite, en vue de 
dStecter vine singularity d 1 absorption et/ou de 
diffusion, les photons hyperdif fusif s non absorbes par 
le milieu sur un large domaine spectral autorisant 

25 1' utilisation des operateurs derives, ces photons 
hyperdif fusif s ou rescap€s Stant €mis par le milieu 
lors de 1 'gclairement de celui-ci. 

De preference, on utilise conjointement le 
temps de vol de 1' impulsion lumineuse et des donnees 

30 spectrales pour Stablir une carte d 1 identity du milieu 
non limpide et les derivees partielles associies. 
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L' exploitation judicieuse de la nature statistique de 
cette carte d'idendite permet de qualifier les echelles 
et les modes d'homogeneisation du milieu sonde vis-&- 
vis des teneurs en absorbeurs et diffuseurs. 
5 On peut eclairer le milieu non limpide par 

une unique impulsion lumineuse S large spectre, 

Selon un mode de mise en ceuvre particulier 
du precede objet de 1' invention, on exploite, en vue de 
qualifier le degrS d'homogSngisation du milieu sondS ou 

10 de detecter une singularity d ! absorption et/ou de 
diffusion, les d£rivees partielles liges aux variations 
en temps de vol et en spectre des photons diffusifs et 
ergodiques non absorbes sur un large domaine spectral 
par le milieu, ces photons diffusifs et ergodiques 

15 etant emis par le milieu lors de 1 1 6clairement de 
celui-ci. 

On peut utiliser conjointement les derivees 
partielles liees au temps de vol de l 1 impulsion 
lumineuse et aux donn£es spectrales pour gtablir une 
20 carte d'identite spectro-temporelle du milieu non 
limpide. 

On peut Eclairer le collecteur sans sonder 
le milieu et en meme temps le milieu non limpide par 
une ou plusieurs impulsions lumineuses autorisant une 

25 imagerie spectro-temporelle derivable par rapport k la 
longueur d'onde et au temps, & deux pics simulant un 
double faisceau. 

La pr§sente invention concerne aussi un 
dispositif de spectrophotom£trie dif fgrentielle, pour 

30 analyser un milieu non limpide, ce dispositif Stant 
caracteris§ en ce qu'il comprend : 
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- une source lumineuse impulsionnelle pour 
eclairer le milieu non limpide et permettant 
1' utilisation ulterieure de la deriv^e partielle, en 
fonction de la longueur d'onde, d'une image spectrale 

5 et temporelle de transmission acquise & partir du 
milieu ainsi SclairS, 

- des moyens d ! acquisition, & partir du 
milieu ainsi eclairg, d ! une image spectrale et 
temporelle de transmission, avec ou sans balayage 

10 spectral, en mode de comptage, et avec un pas en 
longueur d'onde et un pas en temps permettant 
1' utilisation ulterieure des dSrivSes partielles de 
1' image en fonction de la longueur d'onde et du temps 
de vol de 1 ' impulsion lumineuse, et 

15 - des moyens de traitement de cette image, 

consid£r£e comme ■ le moment d'ordre zero, et de ses 
d£riv§es partielles en fonction de la longueur d'onde 
et du temps, pour acqu^rir des informations sur ce 
milieu. 

20 Selon un premier mode de realisation 

particulier du dispositif objet de 1' invention, la 
source lumineuse impulsionnelle comprend des moyens de 
generation non linSaire d ! impulsions lumineuses 
femtosecondes ou picosecondes ou nanosecondes 

25 permettant 1' utilisation de la derivge partielle en 
fonction de la longueur d'onde. 

Selon un deuxieme mode de realisation 
particulier du dispositif objet de 1' invention, la 
source lumineuse impulsionnelle comprend des moyens de 

30 generation et d' amplification de continuum femtoseconde 
ou picoseconde ou nanoseconde d'un continuum, le 
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continuum permettant directement (sans balayage) 
1' utilisation de la dSrivee partielle en fonction de la 
longueur d'onde. 

Les moyens d 1 acquisition peuvent comprendre 
5 une camera & balayage de fente en mode comptage ou tine 
camera k balayage de fente de comptage de photo- 
electron unique, en mode d' operation monocoup ou en 
mode d' operation synchro-scan ou en mode analogique. 

10 On propose en outre un precede de 

spectrophotometrie conforme au proc£d§ objet de 
1' invention et permettant 1 1 gtablissement d'une 
veritable carte d' identity optique du volume sond£, 
c'est-&-dire d'une fiddle signature de la nature 

15 statist ique du milieu vis -S- vis des teneurs plus ou 
moins homogenes des diffuseurs et des absorbeurs, cette 
carte d' identity optique etant sous la forme d'une ou 
plusieurs images spectro-temporelles permettant d' avoir 
accSs simultanSment aux distributions temporelles pour 

20 xine fenetre spectrale donnie, aux distributions 
spectrales pour un intervalle de temps donne, aux 
derivees partielles de ces deux distributions ainsi 
qu'S leurs intigrales. 

On propose aussi un proc§d6 de 

25 spectrophotometrie et tomographie conforme au proc#d£ 
objet de 1' invention, dans lequel on utilise 
conjointement les dSrivSes partielles, par rapport a la 
longueur d'onde, au temps de vol et k l'espace, des 
images spectro-temporelles et des modes monopoints avec 

30 balayage spatial ou des modes multipoints commutes des 
zones d' injection et/ou des zones de collection de la 
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lumidre, 1' utilisation des derivees partielles de type 
espace etant alors possible et permettant de traiter 
certains cas de densite non constante dans le volume 
sondS, ce proc§d£ de tomographie avec des 
5 dif ferentielles d' image spectro-temporelle permettant, 
d'une part, un repgrage d'une singularity de la 
concentration d'absorbeurs et/ou de diffuseurs et, 
d' autre part, une identification molSculaire de ces 
absorbeurs, les injections et collections de la lumiere 
10 pouvant etre r§alis£es soit £. la surface soit au sein 
du volume. 

On propose egalement un procidS de 
spectrophotometrie dif f erentielle conforme au procede 
objet de 1' invention, dans lequel on utilise 

15 con j ointement les dSrivees partielles, par rapport a la 
longueur d'onde et au temps, de l'imagerie spectro- 
temporelle de transmittance et le mode de comptage par 
binarisation de 1' image puis detection d'une zone de 
pixel attribuable a un photo -Electron et reduction de 

20 cette zone k un seul pixel allumS, ou a Tine 6chelle 
subpixel, afin d'augmenter la dynamique et de qualifier 
la mesure en photo- electron unique. 

Selon un mode de mise en oeuvre particulier 
du procSdS objet de 1' invention, on exploite, en vue de 

25 mesurer de tres faibles variations d'un absorbeur, les 
d^rivges partielles en temps et en longueur d'onde de 
la queue temporelle liee aux photons ergodiques 
diffuses par le milieu lors de 1 ' Sclairement de celui- 
ci. 

30 On peut exploiter, en vue de mesurer ou 

d'homog£n€iser des teneurs d'absorbeurs et/ou 
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diffuseurs, l'imagerie spectro-temporelle et les 
opgrateurs d/dt, tf/dt 2 , d/dX, ?/dX 2 , SP/dtdX jusqu'aux 
ordres supSrieurs, ces operateurs s'appliquant & 
1' image spectro-temporelle. 
5 Ainsi, l'imagerie spectro-temporelle 

utilis6e dans la pr^sente invention permet d' avoir 
accSs k des operateurs du type cP/dt m dX n ~ m avec m<n. 

Dans la prSsente invention, on peut 
utiliser con j ointement les fluctuations des images 

10 spectro-temporelles de transmittance et les operateurs 
associ^s (pour les fluctuations sur un temps 
macroscopique par exemple i l'Schelle de la seconde, 
l'operateur d/dt^croacopiqae et les derivees nidmes) afin 
d'effectuer des mesures dynamiques d'opacim^trie, de 

15 colorimStrie et de granulometrie . 

BREVE DESCRIPTION DES DESSINS 

La prSsente invention sera mieux comprise a 
la lecture de la description d'exemples de realisation 

20 donnas ci-apres, a titre purement indicatif et 
nullement limitatif, en faisant reference aux dessins 
annexes, sur lesquels : 

- la figure 1 est une vue sch§matique d'un 
dispositif conforme a l 1 invention, 

25 - la figure 2 illustre schematiquement 

l 1 image spectro-temporelle de transmission en 
enregistrant 1' impulsion blanche avant les milieux a 
sonder et en m§me temps 1' impulsion transmise par ces 
milieux avec le dispositif de la figure 1, 
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- la figure 3 illustre schematiquement 
l 1 image spectro-temporelle de transmission fournie par 
les moyens d'affichage que comporte le dispositif de la 
figure 1, 

5 - la figure 4 illustre schematiquement des 

exemples de la nature statistique des milieux a deux 
phases que l f on peut etudier avec I 1 invention, et 

- la figure 5 est une vue schgmatique d'un 
exemple du dispositif de l 1 invention. 

10 

EXPOSE DETAILLE DE MODES DE REALISATION PARTICULIERS 



On explique d'abord le principe general de 
l 1 invention qui concerne en particulier 1 1 exploitation 

15 de photons « ergodiques > ou « rescapes » pour la 
chimie analytique. 

Dans un milieu limpide, le passage de la 
mesure d'une transmittance optique a une concentration 
(ou & plusieurs concentrations de divers colorants) est 

20 obtenu par la classique loi de Beer-Lambert, qui est 
utilisge largement en spectrophotomitrie . 

L' invention consiste notamment a etendre 
cette approche a des milieux non transparents, par 
exemple des boues industrielles, des milieux contenant 

25 des suspensions, du lait, du fromage, de la viande, des 
tissus biologiques, et plus g£n€ralement a des milieux 
turbides c'est-S-dire des milieux diffusants et 
absorbants . 

L 1 invention peut §tre consid6r£e comme une 
30 generalisation de la spectrophotom§trie classique 
continue; elle peut §tre presentee comme vine technique 
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de colorimetrie et d' opacimetrie des milieux complexes. 

Elle permet aussi de qualifier le degre 
d'homogeneite d'un milieu sondS et de discerner quand 
les grandeurs de concentration moyenne d'absorbeurs 
5 et/ou de diffuseurs ont un sens pour le volume 
d'intSrit (volume consider^ du milieu). 

Avec 1' utilisation des derivees partielles, 
1' invention exploite en particulier les photons 
hyperdif fusif s non absorb^s sur un large domaine 

10 spectral, allant de l'ultraviolet k 1' infrarouge, et 
sur un large domaine de temps de vol, pour detecter une 
singularity d' absorption et/ou de diffusion dans un 
volume constituS de mat6riaux par exemple granulaires, 
poreux ou fractures. 

15 Les Svenements rares de comptage de tels 

photons, const ituant les photons ergodiques, portent 
une information importante de valeur moyenne volumique 
dont le sens statist ique peut notamment §tre compart 
aux taux de survie. II faut noter que l'on mesure les 

20 propri£t6s d' absorption par des photons non absorbes et 
done survivants. 

Les d6riv€es partielles par rapport a la 
longueur d'onde permettent de separer les caracteres de 
monotonie diffSrents entre absorption et diffusion. Les 

25 derivees partielles par rapport au temps de vol 
determinent des critdres object if s et sensibles sur la 
nature du couplage entre les modes diffusifs lies k 
1' organisation macroscopique et m£soscopique du milieu 
et sur les absorptions mol€culaires . Les dSrivSes 

30 partielles croisees permettent de lever les 
incertitudes de conjonction entre les di verses 



WO 2004/013617 



PCT/FR2003/002411 



propagations et absorptions. 

A la fin de la presente description, on 
donne un exemple de m£thode d' utilisation des dirivees 
partielles de 1' image spectro-temporelle. 
5 Si les modes de diffusion restent 

constants, l 1 invention permet de detecter de tres 
faibles variations d' absorption du fait de 
1' augmentation des trajets et de la probability 
d'aborption qui correspondent aux photons consideres. 

10 Pour une distribution constante d'absorbeurs (si le 
milieu est isotrope & de multiples 6chelles, et pour 
une meme density moyenne) , 1' invention permet de 
detecter de faibles variations de diffusion. 

Du point de vue technique, I 1 invention 

15 repose sur l'obtention et l 1 analyse d'une imagerie 
spectro-temporelle de transmission sans balayage 
spectral en vue de 1' utilisation des operateurs de type 
d6rivee partielle. Cette imagerie spectro-temporelle 
differentielle est obtenue par le couplage d'une source 

20 lumineuse impulsionnelle S large spectre, allant de 
1' ultraviolet a 1' inf rarouge, et d'une camera & 
balayage de fente ("streak camera") en mode de comptage 
ou en mode anal ogi que. 

La rSponse impulsionnelle d'un volume 

25 d'int^ret donne accSs k la fonction de transfert 
optique sous la forme d'une image dans un repere dont 
le premier axe est une fen^tre spectrale et dont le 
deuxiSme axe est une f enStre temporelle qui correspond 
H une deflexion temporelle. 

30 L # utilisation conjointe du temps de vol de 

la lumiere et des donn^es spectrales et de leurs 
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deriv6es autorise 1 'etablissement d'une veritable 
"carte d'identite optique" du volume sonde, c'est-a- 
dire d ! une veritable signature de la nature statistique 
du milieu vis-a-vis des teneurs plus ou moins homogSnes 
5 des diffuseurs et des absorbeurs. 

De manidre imagge, en assimilant 
grossi£rement les photons a des « marcheurs », et afin 
de souligner 1 ' importance des aspects statistiques, on 
peut considerer le proced£ optique propose comme suit. 

10 A un instant donn£ des milliards de 

marcheurs sortent d'un stade en un temps trSs bref. 
Chaque marcheur a un maillot de couleur pouvant aller 
du bleu au rouge suivant sa taille, les bleus 6tant les 
plus grands. Chaque marcheur va a une vitesse 

15 « moyenne* de 1'ordre de 0,214mm/ps, les bleus etant 
ISgerement plus lents en moyenne que les rouges. 

Le dispositif de comptage est & une 
certaine distance connue de ce stade. Entre ce point 
d' observation et le stade, se situent les milieux a 

20 etudier, qui, d'une maniSre trds imagee et quasiment en 
deux dimensions, peuvent §tre une ville complexe. 

L'observateur compte le nombre de marcheurs 
suivant leurs couleurs et il cumule ses comptages pour 
une fen§tre temporelle par exemple toutes les 2ps sur 

25 un intervalle de . temps de par exemple 960ps 
(2ps/pixel*480pixels) . Cette s£rie de 480 points est 
directement utilisSe pour approximer les d§riv6es 
partielles en fonction du temps de transport. 

Certains de ces marcheurs sortent 

30 directement du stade i l'endroit oil le dispositif de 
comptage est positionnS. Ce sont des photons dits 
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serpentiles ou balistiques. lis ont ete beaucoup 
etudi^s et ont fait l'objet de nombreux brevets car ils 
autoriseraient en particulier une iraagerie par 
projection comme en radiographic X. Les photons 
5 « ergodiques x> - sont les photons qui ont les temps de 
contact les plus longs avec les milieux a etudier. 
Entre ces deux cas extrimes, de nombreux modes de 
transport sont possibles. Et les derivees partielles 
diagnostiquent la nature de ces modes. Ce n'est pas 

10 seulement la distribution temporelle mais ce sont 
surtout les variations de distribution temporelle de 
ces marcheurs entre les serpentiles et les ergodiques 
qui apportent une information sur la nature statistique 
de 1' exploration des milieux. Les variations de la 

15 distribution spectrale de ces marcheurs renseignent 
surtout sur la nature statistique de la « mort » par 
absorption de ces marcheurs. Finalement, c'est en 
enregistrant les deux distributions simultanement et en 
determinant leurs d§rivees partielles que l'on obtient 

20 la carte d' identite des milieux sondes. 

Dans les milieux limpides, entre le stade 
et le lieu de comptage, on peut reduire le probldme du 
volume sondS a une longueur, d'oti un espace de travail 
dont la dimension est Sgale & 1. 

25 Dans les milieux non limpides et non 

forcSment homog£n6i sables, le probldme demeure de 
nature volumique et sa dimension « pseudo-fractale » 
(c'est-i-dire la dimension de 1' espace de travail) est 
comprise entre 1 et 4 et constitue une inconnue. 

30 L'exemple imag6 de la ville est en fait un 

cas oil la dimension de I 7 espace de travail est comprise 
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entre 1 et 3 si les lieux de transports sont tous sur 
une meme surface. 



simple d'un treillis, le mathematicien G. P61ya a 
5 montre en 1921 que la probability pour que le marcheur 
alSatoire revienne en son point d'origine est egale a 1 
(au bout d'un temps infini) . Ici, 1' ensemble des points 
• d'un treillis est sans aucune absorption. Pour les 
dimensions supSrieures & 2 ce n'est plus le cas. On 

10 obtient les c€l£bres constantes de Polya [voir Glasser, 
M. L. et Zucker, I. J. "Extended Watson Integrals for 
the Cubic Lattices." Proc. Nat. Acad. Sci. U.S.A. 74, 
1800-1801, 1977]. Pour les dimensions 3 et 4, les 
probability valent respectivement 0,340537 et 0,1932. 

15 Ainsi, meme avec un mode de diffusion 

simple (marcheur al6atoire) sur un treillis de 
dimension 3 et m§me avec un temps d' integration trSs 
grand, les photons ergodiques ne peuvent sonder 
totalement le volume. Si le milieu tend vers la 

20 dimension 2, les photons ergodiques peuvent tendre vers 
un sondage quasi -complet pour un milieu tr§s tres 
faiblement absorbant. Malheureusement, les modes de 
structuration des diffuseurs et des absorbeurs sont 
loin de ces cas simples. L' utilisation de la 

25 spectomStrie dif fgrentielle en longueur d'onde et en 
temps de vol, avec des d6placements des points sources 
et/ou des points de detection devrait permettre de 
traiter 1' ensemble des cas reels, qui sont loin d'un 
simple treillis de dimension entiere. 

30 Un dispositif correspondant & ce qui 

prScSde est sch£matiquement illustre sur la figure 1 oti. 



Pour les dimensions 1 et 2 et dans le cas 
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l'on voit un volume d'int£r§t V, une source lumineuse 
impulsionnelle a. large spectre S, un point d 1 impact P 
d'une impulsion emise par cette source sur le volume V, 
des moyens 2 de detection de la lumiere ayant traverse 
le volume V, des moyens electroniques 4 de traitement 
des signaux fournis par les moyens 2 qui sont 6galernent 
appeles "collecteur" (et qui peuvent etre, par exemple, 
une fibre optique, une lentille ou une micro-cam§ra) et 
des moyens 6 d'affichage dont sont munis les moyens 4. 

La figure 2 illustre schSmatiquement 
1 1 image spectro-temporelle de transmission en 
enregistrant 1' impulsion blanche avant les milieux a 
sonder et en meme temps 1' impulsion transmise par ces 
milieux. 

On y voit la fonction instrumentale 
(suivant 1' image des marcheurs, la convolution de la 
sortie du stade et de la detection sans passer par la 
ville, dans l'air par exemple) avant le milieu (I) et 
aprSs le milieu (II) . On a portS en abscisse la 
longueur d'onde X de la lumiere (en nanometres) et en 
ordonnee le temps t (typiquement en picosecondes par 
pixel) . 

Le dispositif et le proc£dS objet de 
1 ' invention permettent 1' enregistrement de la fonction 
instrumentale et de la rgponse optique au travers des 
milieux sondes. Get enregistrement est alors realist en 
m€me temps, simulant un double faiseau par porte 
temporelle adapt 6e, assez espac§e pour eviter une 
accumulation (« pile -up ») . 

La figure 3 illustre schSmatiquement 
I 1 image spectro-temporelle de transmission fournie par 
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les moyens d'affichage 6, avec 1' enregistrement de 
seulement la fonction instrumentale ou seulement la 
r£ponse optique au travers des milieux sondes. La 
longueur d'onde X (en nanometres) est portee en 
5 abscisse et le temps t (typiquement en picosecondes par 
pixel) en ordonnee. On prScise que I 1 image a 
correspond, dans l'exemple represents, au comptage de 
photo-61ectrons en mode analogique. 

Soit I 0 l 1 intensity de la lumi§re incidente 

10 qui va traverser le volume V, suivant la notation 
classique utilisee en spectrophometrie de Beer-Latnbert 
des milieux limpides. Dans cette spectromStrie non 
impulsionnelle, il est tres courant d' avoir des 
dispositifs avec deux faisceaux, l'un passant par la 

15 cuve coloree et 1' autre permettant de mesurer ce I 0 . 

L' invention permet aussi de realises une 
imagerie spectro-temporelle deux coups" grace au 

caractdre impulsionnel : on realise la mesure de I Q ou 
de la fonction instrumentale, c'est-a-dire de la 

20 convolution de la sortie du stade et de la detection 
sans les milieux k Studier, vis-&-vis des analyses du 
temps de vol des photons. On precise que la mesure de 
I 0 peut §tre r^alisSe gr3.ce a une fuite de lumiere, 
avant la propagation dans le volume V, ou de fagon 

25 sSquentielle. 

On precise qu'en mode inf£rieur & 1000 
f emtosecondes , 1'Snergie envoySe est contenue dans une 
impulsion de longueur Squivalente inferieure a 0,3mm 
dans le vide. Mais si 1'objet a analyser a des 

30 grandeurs caracteristiques de l'ordre d f une centaine de 
metres, on peut utiliser des impulsions nanosecondes 
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comme celles qui sont fournies par les lampes-f lash & 
hydrogene, et des temps de deflexion d'une centaine de 
nanosecondes . 

Cependant, pour des raisons de probability 
5 d' absorption et surtout de niveau de comptage, 
l 1 invention est surtout utile pour sonder des volumes 
inf^rieurs a une dizaine de m 3 jusqu'S. des volumes de 
1'ordre du mm 3 . 

La limite inferieure du volume que l'on 

10 peut sonder ne depend que des performances du systdme 
de deflexion de la camera a balayage de fente, sur 
laquelle on reviendra par la suite, car des sources & 
large spectre (et qui correspondent a des sources 
ultra-breves de quelques fs) sont dSj& connues. 

15 On decrit maintenant un exemple du procSdg 

objet de 1' invention. 

On commence par injecter, soit a la surface 
soit dans le volume de 1' objet & mesurer, une impulsion 
ultracourte, k tr£s large spectre dans la fenetre 

20 constitute par les domaines ultraviolet, visible et 
proche-infrarouge. Les exemples les plus connus de ce 
type d' impulsion sont les impulsions f emtosecondes 
ultracourtes a large spectre en sortie d' oscillateur , 
les impulsions " resultant de la g6n6ration d'un 

25 continuum femtoseconde ou picoseconde et/ou de 
1' amplification d'un continuum. 

L' injection et la collection de la lumiere 
peuvent etre r6alisees soit en mode monopoint soit en 
mode multipoint avec ou sans dSplacement automatique. 

30 L'utilisation d'une pile d' images spectro- 

temporelles avec ce genre de deplacement permet une 
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tomographic optique c'est-&-dire (1) un reperage d'une 
singularite de la concentration d'absorbeurs et (2) une 
identification moleculaire de ces absorbeurs. 

La collection peut Stre r^alisee soit k la 
5 surface soit au sein du volume. 

La lumiere collect£e est analysee par un 
§16ment diffractant (spectrographe a r£seau 
(« grating ») ou a prisme ou autre appareil de 
diffraction optique) permettant 1' utilisation de la 

10 dSrivee en fonction de la longueur d'onde et couple k 
une camera a balayage de fente ou a tout appareil 
permettant de deflechir temporellement la lumiere d ! une 
fenetre spectrale autorisant 1' utilisation de la 
derivee en fonction du temps de vol. 

15 On obtient ainsi une image spectro- 

temporelle de *1' arc-en-ciel" deflechi temporellement, 
qui est derivable. Diff brents modes d' analyse 
mathSmatique peuvent §tre mis en oeuvre pour resoudre 
les problemes inverses qui visent a quantifier des 

20 grandeurs homogen6isables du milieu traverse. 

Des aspects essentiels de 1' invention sont 
1' utilisation de cette image et de ses derivees 
partielles en tant que moyen qualifiant et quantifiant 
1' ensemble des propriStSs optiques d'un objet et 

25 l'arriv6e k divers degres de formation de moyenne 
volumique ou la decision de ne pas former une moyenne 
volumique . 

Rappelons que, pour une longueur d'onde 
donnee, on dSfinit g§n§ralement trois types de photons 
30 suivant la distribution temporelle: 

a) les photons balistiques qui sont les 
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premiers arrives et n 1 ont pas connu un €vdnement de 
diffusion, 

b) les photons serpentiles qui arrivent 
juste apres les photons balistiques et restent 

5 utilisables a des fins d' imager ie dans des milieux 
turbides, et qui malgre quelques SvSnements de 
diffusion gardent une information de leur direction 
initiale, et 

c) les photons diffusifs qui arrivent bien 
10 apr§s les precedents. 

Dans la presente invention, on considSre 
aussi une quatrieme classe, a savoir la classe des 
photons hyper-diffusifs ou ergodiques. Malgr§ leur long 
temps de contact avec le milieu (ou long temps de vol) , 

15 ceux qui sont mesures ont &t£ non absorbSs comme les 3 
autres classes. 

Au del§. de ces 3 ou 4 classes, 1' invention 
se concentre sur les dSrivSes partielles en fonction de 
la longueur d'onde et du temps de vol et sur le mode 

20 comptage. 

La mesure de cette classe de photo- 
electrons uniques (SPE) hyperdif fusif s correspond k des 
gvenements trds rares. Ces photons portent une 
information considerable et souvent decisive. 

25 Malheureusement , il faut relativiser leur importance £ 
cause de leur faible taux de comptage que permettent de 
compenser les modes de comptage. 

De plus, suivant 1' image du marcheur, il 
reste difficile de savoir si ces photons ergodiques ont 

30 sonde 1 # ensemble du volume d'int£ret ou si ils sont 
rest£s piiges dans une singularity hyperdif fusante (cas 
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tr£s eloigne du treillis ideal de Polya) . Par contre, 
si le milieu se rapproche du treillis isotrope 
pr^sentant une isoprobabilitS, la moindre fluctuation 
de la teneur des absorbeurs leur est souvent fatale. 

5 Mais, dans 1 ' ensemble, la distribution 

spectrale, les variations de cette distribution 
(dSriv§e premi&re) et les d6riv6es ni&mes et les 
d£rivees crois£es de 1' ensemble des photons et plus 
particulidrement des photons hyper-dif fusif s permet de 

0 lever ces types d' ambiguity. La dynamique de la mesure 
est done cruciale. Le passage en SPE acroit 
considerablement cette dynamique. Outre les axes 
longueur d'onde et deflexion temporelle, l'axe z en 
unites SPE est tres important. L'imagerie spectro- 

5 temporelle en mode comptage pr£sente un tres grand 
potentiel car elle atteint les limites de ces 
techniques optiques vis-a-vis de la dynamique de mesure 
surtout pour compter et accumuler ces photons 
ergodiques . 

0 L' image des marcheurs est dangereuse. En 

fait, la propagation de la lumiere dans un milieu dense 
est consid§r§e comme une perturbation d'un champ 
d' induction electrique et d'un champ d' induction 
magn6tique, tous deux oscillant avec une pSriode de 1 a 

5 3 femtosecondes ■ (par exemple pour la fenStre 
300-900nm) . La propagation de la lumiere est alors 
r^alisee par des dipdles, de tailles moleculaires, qui 
rSinduisent un champ Slectromagnetique de manidre 
amort ie. La structure du milieu dense a des Schelles 

0 supramolSculaires induit dif fSrents regimes de 
diffusion/diffraction (Rayleigh, Mie, Fraunhofer) . 
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Dans un milieu limpide, des phenomenes de 
diffusion ont lieu mais n'entralnent qu'une variation 
de direction (loi de la refraction). L'6cart-type de 
cette direction depend des fluctuations locales de la 
5 densite et de 1' amort is sement des dipoles c'est-&-dire 
de la distribution de l'indice de refraction n qui peut 
§tre souvent approximee par une distribution de Dirac. 
Par exemple, pour une temperature homogSne fixge et 
pour une longueur d'onde donnSe, n vaut 1,4. 

10 On peut considSrer tous les milieux 

diffusants comme un milieu polyphasique non periodique 
ou diverses echelles plus ou moins proches de la 
longueur d'onde interviennent . Dans ce type de milieu, 
les ondes secondaires c'est-a-dire remises par les 

15 dipdles oscillants n' interf &rent plus de manidre 
constructive pour privilSgier une direction. 

Contrairement aux ondes quasi -directes 
(souvent nomm6es photons balistiques ou serpentiles) , 
les ondes diffusives ont beaucoup moins retenu 

20 1' attention. Le cas des « photons ergodiques » fait 
qu'ils peuvent §tre consid^res comme des ondes 
r^siduelles avec un tr£s fort dephasage ou un §norme 
retard. Ces photons ergodiques mesures ont St§ de fait 
r§£mis par de nombreux dip61es oscillants non 

25 absorbants. 

Ces photons hyperdif fusif s ont 6te souvent 
consid§rfis comme 6tant les moins porteurs 
d' information. 

S*il est vrai qu'ils peuvent etre 
30 considSres comme pauvres du point de vue de l'optique 
geometrique et des techniques d'imagerie, ils sont 



WO 2004/013617 



PCT/FR2003/002411 



neanmoins riches d'une information de moyenne volumique 
au sens des marcheurs ou aussi au sens des ondes 
electromagnetiques. Par exemple, ces photons 
hyperdif fusif s sont de loin les plus sensibles & la 
5 moindre singularity absorbante et plus globalement & la 
partie imaginaire des indices de refraction. 

L'auteur de la presente invention a trouve 
que ces photons hyperdif fusif s sont susceptibles de 
former un excellent moyen pour sonder un volume 

10 complexe et portent beaucoup d' informations sur ce 
volume explore raais qu'il faut des degres de liberty 
suppl£mentaires pour decoupler 1" absorption et les 
diff brents regimes d' induct ion/ diffusion que permet la 
spectromgtrie dif f Srentielle en spectre et en temps de 

15 vol, Le mode de comptage met aussi en valeur cette 
analyse . 

La pertinence de la reponse impulsionnelle 
a large spectre, et non k quelques longueurs d'onde 
discretes, repose en partie sur le fait que 
20 -les diffusions sont le plus souvent 

monotones en fonction de la longueur d'onde et 

-les absorptions sont le plus souvent non 
monotones (existence de pic ou de trouj , comme cela est 
bien connu en spectrophotomStrie dif f erentielle et en 
25 chimiomStrie optique. 

Dans les classes de la distribution 
temporelle des photo- Electrons mesurSs, celles qui 
portent le plus d' information vis-a-vis du volume sonde 
sont les derniSres classes, c'est-a-dire celles des 
30 temps longs. 

De plus, la probability de non-absorption 




WO 2004/013617 PCTYFR2003/002411 

25 

au bout d'un "trajet" de L mm varie comme exp (-L/La) 
oil La (en mm) est la longueur d' absorption, ce 
qu'exprime la loi de Beer-Lambert. Plus le trajet 
augmente, plus la probability de non absorption diminue 
5 de mani£re exponentielle.Dans le cas de 1'imagerie 
proposee dans la presente invention, on peut aussi 
parler d'imagerie logarithmique . La spectromStrie 
dif f6rentielle permet vine utilisation judicieuse de ce 
caractere logarithmique . 

10 . En outre, le passage du "temps 

d' interaction" i la "longueur d' interaction" conduit & 
l'ecriture de la vitesse moyenne. 

Le problSme geom£trique est en fait plus 
complique. Ainsi, il n'est pas aisS de separer un 

15 photon hyperdiffusif , qui est reste longtemps dans un 
petit volume en hypersondant ce dernier, de celui qui a 
sond§ un volume plus grand et qui n'est pas passe deux 
fois sur le mime site (dSnommS trajet auto-excluant) . 
On peut seulement af firmer que I'on connait les points 

20 d'entrSe et de sortie et son temps d' interaction. 

Remonter a un trajet d' interaction 
necessite des hypotheses ou des connaissances sur la 
structuration du milieu. 

On precise que l f invention ne se limite pas 

25 a 1 'utilisation conjointe du spectre de transmission et 
de la distribution du temps de vol et de leurs d£rivees 
partielles, Elle concerne aussi l'utilisation d'un 
dSplacement des points d' impact et/ou de collection 
afin de rSaliser une tomographie (mettant en oeuvre des 

30 deriv§es partielles de type espace) . Avec ces degrgs de 
liberty supplimentaires, les distributions spectrale et 
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temporelle permettent de mieux qualifier et situer dans 
un volume des singularites de teneur d'absorbeurs et/ou 
de diffuseurs. 

Plus generalement 1 ' imagerie spectro- 
5 temporelle permet d 1 avoir accds a 1'operateur <f/dtd% et 
aux ordres superieurs, avec la fonction peigne 
*d'§chantillonnage (par exemple 480 pixels pour 1,1ns et 
640 pixels pour I80nm, 6tant entendu que cela peut 
changer suivant le polychromateur et la camera utilisee 

10 au bout de la chalne d' acquisition) . 

En spectrophotomStrie dif f€rentielle 
classique, on n f a accSs qu'aux opirateurs d n /dX n , et, 
avec les techniques de temps de vol a quelques 
longueurs d'onde discretes, aux opirateurs d n /dt n . En 

15 chemometrie classique, on ajoute les opSrateurs dP/dCf, 
ou Ci est la concentration du i-idme absorbeur, et les 
operateurs croisSs du type d 2 /dCid%. L' imagerie propos6e 
dans la pr^sente invention permet d' avoir accds & de 
nouveaux operateurs. 

20 DiffSrents modes d' analyse mathematique 

peuvent, Stre mis en oeuvre pour r^soudre les probldmes 
inverses qui , visent a quantifier des grandeurs 
homogenSi sables du milieu traverse. 

On poura consulter le document suivant : 

25 Kienle A., Patterson M. S., (1997) improved solutions 
of the steady state and the time-resolved diffusion 
equations for reflectance from a semi -infinite turbid 
medium, J. Opt. Soc. Am. A 14 246-254 

-pour une ggometrie semi-infinie 

30 -pour un espace de travail dont la 

dimension est posee ggale 3. 3, 
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-pour un point d' impact infiniment bref et 
d'un diamStre infiniment petit sur la surface, 

-pour un point de collection sur la surface 
assez loin du point d' impact, et 
5 - pour un milieu assez homogine od 

1' approximation de la diffusion de la lumidre est 
valide et ou les teneurs des absorbeurs et des 
diffuseurs sont homogenSisables et rgductibles a deux 
parametres (vu en mm" 1 et ]i 3 ' en mm" 1 ) . 

10 A la fin de la prSsente description, on 

donne quelques exemples de traitement et de l'int6ret 
des nouveaux operateurs. 

De nombreuses approches mathematiques ont 
6t& proposees dans l'art anterieur pour r^soudre- le- 

15 probleme inverse des milieux turbides, aussi bien du 
point de vue de la mesure des moyennes des 
concentrations des absorbeurs ' et des diffuseurs que du 
point de vue de la moyenne des g€om6tries des 
diffuseurs (granulometrie) . 

20 L'imagerie spectro-temporelle telle que 

dScrite permet d' avoir acc£s simultanSment aux 
distributions temporelles et spectrales et aux derivSes 
partielles de celles-ci ainsi qu'aux int£grales de 
celles-ci. Ces moyens permettent de qualifier la nature 

25 statistique des milieux traverses (optiquement 
« moyennisables » ou non) . Si des connaissances 
gSometriques ou de diverses teneurs locales des 
absorbeurs et diffuseurs sont ajoutSes, le probleme 
math£matique est simplifiS et peut permettre, suivant 

30 les cas, de mesurer plus pr6cis6ment les teneurs des 
absorbeurs et des diffuseurs. 
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Quels que soient les choix des mSthodes 
d' analyse mathematique (methodes de type Monte Carlo, 
equations radiatives, ou Equations maxwelliennes par 
exemple) , l'imagerie spectro-temporelle proposee dans 
5 la presente invention procure une carte d'identite 
optique et empirique des milieux sondes. L'6tude 
statistique de ses variations par 1' intermediaire des 
derivees partielles perraet de qualifier directement les 
variability des milieux explores. 

10 En outre, 1' exploitation des photons 

hyperdif fusif s permet le passage de 1' image spectro- 
temporelle a des informations quant if iees au sujet du 
volume Studi€. Ces formations sont soit de nature 
gSometrique (tomographic, detection de singularity) 

15 soit de nature chimique (chemometrie, quantification de 
concentration) . 

L' invention se situe ainsi k deux niveaux: 

1) €tablissement de la carte d' identity 
optique (images spectro-temporelles et leurs derivSes 

20 partielles en longueur d'onde et en temps de vol) : on 
se limite k un controle de qualite, 

2) exploitation de la carte d' identity 
optique : on cherche & analyser le parcours des photons 
hyperdif fusif s pour homogSneiser les teneurs des 

25 absorbeurs et/ou des diffuseurs ou decider d'Stablir 
une moyenne. 

Si I'hoiriogeniisation est valide, la 
concentration des absorbeurs et des diffuseurs a un 
sens et est quantifiable, 
30 Revenons de fagon plus dStaillSe sur l'un 

des problemes que r^sout la prSsente invention. 
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On souhaite faire une mesure globale et in 
situ des concentrations ou des variations de 
concentration d'absorbeur et de diffuseur dans un 
milieu non limpide. De plus, 1' invention vise a savoir 
5 quand il est valide de parler de concentration moyenne 
et permet de tester l'hypothdse de 1 'homogen^isation. 

Illustrons briSvement certaines classes de 
problSmes. On peut vouloir: 

-mesurer une concentration de colorant dans 
10 un milieu turbide, par exemple les flavines endogSnes 
dans du lait, ou mesurer des colorants dans une 
suspension, 

-mesurer les concentrations de divers 
colorants, par exemple un colorant bleu, dans du 
15 fromage, 

-dStecter le volume de morceaux de fruits 
rouges (ou d'autres couleurs) dans des pots de yaourt, 
dans une chaine de fabrication, 

-quantifier la composition corporelle par 
20 vine analyse combinee de la repartition de l'eau et des 
graisses, 

-quantifier la quality d'un melange, et 

-d^celer un polluant dans le brouillard. 

Contrairement au cas d'un milieu limpide 
25 contenant un seul colorant, la seule connaissance de la 
ggom6trie des points d' injection/collection de la 
lumiSre, associ^e a la connaissance de 1'intensite 
initiale, ne suffit plus. II faut d'autres donnees. 

Dans le cas des milieux limpides, 1' analyse 
30 des propriites spectrophotom^triques permet de rempnter 
& la concentration suivant la loi de Beer-Lambert. 
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Les appareils connus de spectrophotomStrie 
reposent tous sur cette analyse. 

— Dans le cas de plusieurs colorants dont le 

melange est inconnu, ou dans le cas de plusieurs 
5 colorants dont I'un est inconnu, diverses m£thodes 
chimiomStriques permettent de resoudre certaines 
classes de ces problemes inverses. 

Les mSthodes de chimiomStrie spectrale 
peuvent s'appliquer a de nombreux domaines, en 
particulier en biologie. 

La present e invention propose une 
chimiomStrie spectrale et temporelle en mode 
impulsionnel et dif f Srentiel . 

Dans l 1 invention, on prend en compte les 
objets qui sont diffusants ou turbides et/ou d'une 
complexity qui empSche 1' analyse par l'optique 
gSometrique classique' (loi de Beer-Lambert) . 

Dans ces cas, les trajets lumineux ne sont 
plus assimilables k une distribution approximee par une 
distribution de Dirac. Les trajets de l'onde 
Sleet romagnetique forment une distribution plus ou 
moins complexe qui est le plus souvent une inconnue. 

A ce sujet, on considdre trois types de 

problemes : 

1) vision dans des milieux troubles 
En fait, tout milieu diffusant est un 
syst&me polyphasique avec des dimensions proches de la 
longueur d'onde de la lumidre avec laquelle on Sclaire 
le milieu. 

A ce propos, consid6rons un exemple de six 
milieux a 2 phases. Cet exemple est sch§matiquement 
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illustre par la figure 4 ou les deux phases sont 
respectivement constitu6es par une matrice I et par un 
element II. 

La taille de 1 Pigment est notSe D. 
5 Sur la figure. 4, on a considSr6 six cas 

numerotes de A a F : 

A : element plus transparent et de m§me 
absorption que la matrice, 

B : Aliment plus transparent et de plus 
10 forte absorption que la matrice, 

C : element diffusant comme la matrice et 
de plus forte absorption que la matrice, 

D : element plus diffusant que la matrice 
et de meme absorption que la matrice, 
15 E : Sliment plus diffusant que la matrice 

et de plus faible absorption que la matrice, 

F : 61£ment plus diffusant que la matrice 
et de plus forte absorption que la matrice. 

II faut en plus traiter les differents cas 
20 des frontiSres de ces volumes (conditions aux limites 
absorbantes, reflectives, d'un troisieme type ou 
inconnues) soit entre la matrice et l'exterieur, soit 
entre la matrice et le sous-volume d'interet. 

Si le milieu n'est pas homogeneisable (D 
25 tr§s supSrieur a la longueur d'onde et/ou periodicite 
et/ou periode du volume d'interet trop faible), alors 
le but est de realiser une tomographic, en particulier 
de reperer un ou des sous -volumes dont les propriStes 
optiques sont trds differentes de celles de la matrice. 
30 L' exploitation de l'imagerie spectro-temporelle et de 
ses derivSes partielles necessite alors des 
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displacements entre le point d* impact du faisceau laser 
de mesure et la detection permettant 1' utilisation des 
derivSes partielles de type espace. 

2) spectroscopie dans des milieux troubles 
5 L'imagerie spectro-temporelle permet de 

resoudre les problemes de chimiometrie dans des milieux 
non transparents mais homog£nSi sables au sens 
mathSmatique. Pour la theorie de 1 'homogenSisation en 
particulier dans des milieux mod£lisables par deux 
10 moyennes de manidre periodique, il est recommande de 
consulter le livre suivant: 

N. Bakhvalov et G, Panasenko, 
"Homogenisation: averaging processes in periodic 
media \ 1989. 

15 Ces probldmes de milieux polyphasiques se 

rSvdlent concerner la majorite des cas de mesures sur 
des objets r6els en chimie, en biologie mais aussi dans 
de nombreux secteurs industriels (poudres, aerosols, 
suspensions par exemple) . 

20 Dans ce cas 2), le milieu d'etude peut §tre 

ramene a un ensemble periodique ou comportant un grand 
nombre de cellules (ou de volumes el£mentaires 
reprSsentatif s) du genre des six exemples A I F de la 
figure 4 . 

25 La model isat ion de la geometrie et des 

phenomenes optiques sont l'un des f acteurs-clSs de la 
resolution des problemes a traiter. On n'a consid£r6 
ici que des moddles simples, dits 3. une seule 6chelle, 
mais les diff6rents regimes de 

30 diffusion/diffraction/inductions nScessitent chacun une 
echelle adaptSe. 
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3) deux types de distributions 

La mesure conduit & une image spectro- 

temporelle qui decrit deux types de distributions: 
-une distribution en temps derivable 
5 -une distribution en longueur d'onde 

derivable . 

L' unite est un nombre entier de photo- 
electrons uniques ou en unite arbitraire (pixels 
allumSs) . 

10 Pour les milieux turbides (absorbants et 

tres dif fusants) , l'analyse de ces distributions en 
temps repose souvent sur 1' approximation de la 
diffusion. 

Les classes de distribution spectrale 

15 dependent surtout des spectres d' absorption des 
colorants. Des programmes de calcul en chimie quantique 
permettent d'obtenir les spectres d' absorption de 
molecules assez simples. Sinon, des bases de donnies de 
spectres ultraviolet- visible-proche inf rarouge peuvent 

20 etre utilisees. 

A ce sujet, on se reportera a : 
H.Du et al., "PhotochemCAD: A computer- 
aided design and research tool in photochemistry," 
Photochemistry and Photobiology, 68, 141-142, (1998), 

25 Les problimes oil la distribution spectrale 

est non monotone (existence d'un pic ou trou 
d' absorption ou d'un pic ou trou de transmission dans 
la fen§tre spectrale sond£e) forment la classe la plus 
facile a traiter avec des d£riv6es partielles d'ordre 

30 un et deux. 

Contrairement au cas d'un milieu limpide 
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contenant un ou quelques colorants bien melanges , les 
problSmes que l'on desire resoudre nScessitent, en plus 
de la distribution spectrale, la distribution 
temporelle de 1' impulsion optique et de ses dSriv£es. 
5 Un aspect de 1' invention est la recherche 

de f en§tres d' integration (des distributions spectrale 
et temporelle) permettant d'optimiser le rapport signal 
sur bruit et de faire intervenir la dynamique entre 
comptage de photons diffusifs et de photons ergodiques. 
10 Dans I'art anterieur, les aspects de 

spectrometrie dif f 6rentielle ont ete abordes vis-S-vis 
de la distribution spectrale sans mode impulsionnel . De 
m§me, des etudes sur les distributions temporelles k 
des longueurs d'onde discrdtes ont et£ faites. 
15 La conjonction des apports de la 

spectromitrie different ielle et de 1' analyse de la 
migration des photons const itue un aspect important de 
l 1 invention. 

On donne ci-apr^s un exemple d'un 
20 dispositif conforme a l 1 invention, permettant d'obtenir 
1 1 image spectro- temporelle . 

En ce qui concerne la source 
impulsionnelle, au moins deux choix techniques sont 
possibles : 

25 a) line source d 1 impulsions f emtosecondes 

ultracourtes en elles-memes- (oscillateur ultra-court) , 
et 

b) des moyens de generation de continuum 
femtoseconde ou picoseconde et/ou d' amplification 
30 param£trique d'un continuum monofilament. 

Le detecteur est une camera & balayage de 
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fente. Au moins deux modes de fonctionnement sont 
possibles: 

i) le mode de balayage synchrone ou mode 
synchro -scan 
5 ii) le mode mono-coup. 

La figure 5 il lustre schematiquement un 
exemple correspondant au cas b)-i). 

Le dispositif de la figure 5 comprend 
successivement : 
0 - un laser femtoseconde 8, 

- des moyens 10 de mise en forme de la 
lumiSre fournie par ce laser 8, cette mise en forme 
etant par exemple une simple lentille, la lumidre ainsi 
mise en forme parvenant alors & un volume 11 du milieu 

5 non limpide etudie . pour y engendrer un continuum 
femtoseconde autorisant ulterieurement les derivations 
en fonction de la longueur d'onde, 

- des moyens 12 de filtrage de la lumiSre 
issue du volume 11, 

0 - des moyens 14 de guidage de la lumiere 

ainsi filtree, par exemple une fibre optique, 

- ion polychromateur 16, qui regoit la 
lumidre ainsi guidSe, autorisant ulterieurement les 
derivations en fonction de la longueur d'onde, 

5 - une camera §. balayage de fente 18 qui 

capte la lumiere issue du polychromateur 16 et 
fonctionne en mode de balayage synchrone, autorisant 
. ulterieurement les derivations en fonction de . la 
longueur d'onde et du temps de vol, 

0 - des moyens elect roniques et informatiques 

20 de traitement des signaux eiectriques fournis par la 
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camera 18, ces moyens electroniques et inf ormatiques de 
traitement 20 etant prevus pour assurer la binarisation 
et le comptage des photo-electrons uniques et pour 
fournir 1' image finale et ses derivSes partielles par 

5 rapport au temps de vol et a la longueur d'onde, et 

- des moyens 22 d'affichage des resultats 
obtenus gr&ce aux moyens electroniques et informatiques 
de traitement 20. 

Le temps de vol est dSfini comme le temps 

0 de propagation de 1' impulsion lumineuse entre le point 
d'entrSe dans le volume 11 etudie et le point d' entree 
dans la camera Sl balayage de f ente 18, 

Le dispositif de la figure 5 fournit une 
imagerie sur, par exemple, 640x480 pixels et dont 

5 1' image est d6rivable en longueur d'onde et en temps de 
vol, A titre d' exemple, chaque pixel est cod§ sur 8 
bits puis stocke sur 16 ou 32 bits. Ce nombre entier 
represente le nombre de photo-Slectrons uniques ou de 
pixels allum6s qui sont comptes pendant un certain 

0 temps de mesure dit temps d' integration. 

Le mode analogique est Sgalement utilisable 
mais le mode de comptage en mode "streak camera" est 
plus avantageux, notamment & cause de la grande 
dynamique et du rapport signal/bruit Sieve qu f il permet 

5 d'obtenir. 

II convient de noter que 1' imagerie 
spectro-temporelle de fluorescence connue dans l'etat 
de la technique, par exemple du genre de celle qui a 
6t£ dSveloppee par la sociStS Hamamatsu, donne acces a 
0 des distributions spectrales et temporelles qui sont 
spScif iques de 1' analyse des temps de declin de 
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fluorescence multicomposante. 

L'imagerie spectro-temporelle de 

transmission, qui est mise en ceuvre dans la pr^sente 
invention, donne acces, quant a elle, S des 
5 distributions spectrales et temporelles et a leurs 
dSrivges partielles, qui sont specif iques de 1' analyse 
des modes de propagation/absorption dans des milieux 
non limpides. 

Le probldme de 1' etude de 1 'elargissement 

10 par augmentation de la diffusion de 1' impulsion 
transmise et de 1 1 glargissement par diminution de 
1' absorption de 1' impulsion transmise et en particulier 
du front descendant est aussi de nature statistique. 

Pour rSsoudre les deux classes de problemes 

15 qui se posent (test d'homogeneisation-sens de la 
moyenne d'une concentration, calcul de la concentration 
moyenne et de ses fluctuations), on utilise 1' ensemble 
des informations brutes dont on dispose et surtout de 
1' ensemble des dirivees partielles. 

20 Pour une distance moyenne entre le point 

d 1 impact p (figure 5) de la lumiere emise par le laser 
8 et le d§tecteur, c'est-S-dire la camera 18, ces 
informations brutes consistent en vine image spectro- 
temporelle de transmission derivable. 

25 Contrairement & une image topographique, on 

a alors une grandeur trds diff£rente pour chaque 
dimension, k savoir des longueurs d'onde (en nm) et un 
1:emps de deflexion (en ps) . 

Les distributions sont done trSs 

30 diff6rentes. 

Contrairement k une image spectro- 
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temporelle de fluorescence dont le but est de 
determiner les temps de dSclin de fluorescence de 
divers fluorophores [S. MOTTIN, C. TRAN-MINH, P. 
LAPORTE, R. CESPUGLIO et M. JOUVET, Fiber optic time- 
5 resolved fluorescence sensor for in vitro serotonin 
determination, Applied Spectroscopy, 1993, 47, 590- 
597], l'imagerie spectro-temporelle de transmission 
presente des distributions plus complexes. 

A titre purement indicatif et nullement 
10 limitatif, le laser 8 utilise a les caracteristiques 
suivantes : 

oscillateur Ti:Sa800nm 78MHz 500mW ; chalne 
CPA (Chirp Pulse Amplification) Ti:Sa pompee par un YLF 
10W, permettant d'obtenir 0,7W, IKHz, 150fs k 800nm ; 
15 generation de continuum. 

Avec ces caracteristiques, la cadence est 
de 1000 impulsions par seconde. 

On peut utiliser une camera & balayage 18 
pour laquelle le temps de chaque image est une 
20 integration sur 33ms, ce qui correspond a 33 tirs 
laser. 

Cependant un mode mono-coup et mono-tir 
reste accessible. 

Deux modes de traitement final du comptage 
25 sont possible : 

1) un seuillage et une binarisation (1 
photo-electron donne en moyenne moins de 5 pixels 
allum£s), l f unite etant alors le "count" , et 

2) un seuillage puis une binarisation et 
30 une analyse morphologique visant & rSduire le nombre N p 

de ces pixels allum6s S 1 pixel "central" allumS, 
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1'unitS etant alors vraiment le photo- electron unique. 

Ce deuxieme mode necessite un comptage 

adapts . 

Ces deux modes existent dans 1'imagerie 
5 spectro-temporelle de transmission derivable. 

Pour chaque image (correspondant & une 
integration sur 33ms) , le taux de comptage peut it re 
approximativement calcule: 

Si la distribution est uniforme et si 1'on 
10 se limite a 1% de ratio de comptage pour Sviter que 
deux photo- electrons tombent sur le m§me groupe de 
pixels et soient comptSs pour un seul photo-electron, 
alors on obtient une valeur grossiSre qui est de 
l'ordre de 307200x0, 01«3000 photo-electrons uniques par 
15 image . 

Si la distribution est du type 
" exponent ielle x fonction d'Heaviside" et si 1'on 
desire vine dynamique finale de 1000, avec un ratio de 
comptage [1000-1] sur 415 pixels/480 pixels, on a 
20 alors, sur une deflexion totale de llOOps, une 
distribution De (t) qui est uniforme sur les 640 pixels 
et telle que, pour toutes les longueurs d'onde : 

De (t) =1000xexp (-t/tm) avec tm=1100/8=137ps . 
Si l f on se limite & un ratio de comptage de 
25 1% pour la zone de comptage final [1000-300] , alors on 
a 72 pixelsx640=46080 pixels et 46080x0, 01«460 photo- 
electrons uniques par image. 

Si la distribution spectrale D*(x) n'est pas 
uniforme mais gaussienne et centree sur le 320idme 
30 pixel et si elle est assez large pour avoir une 
dynamique de 1000 sur 52 pixels, alors on peut Scrire 
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D*. (x) =exp ( - ( (x-320) /100) 2 ) ou x est le numero de pixel. 

On pourrait aussi utiliser une distribution 
pour 1' absorption de type lorentzien ou une 
distribution non symStrique comrae la distribution de 
5 Lowry et Hudson. 

A ce sujet on pourra se reporter au 
document suivant : 

C. Sandorfy, "lies spectres electroniques en 
chimie theorique, Paris, Revue d'optique th6orique et 
10 instrumental, 1959. 

De meme, on obtient 

72 pixelsx220 pixels x0,01 «158 photo-electrons uniques 
par image. 

Avec ces deux distributions, le comptage 
15 s§curis6 est tres grossierement de 1'ordre de 160 
photo-electrons uniques par image. Avec 30 images par 
seconde, on arrive & environ 5000 SPE par seconde et en 
moyenne 5 SPE par tir laser & la cadence de 1kHz. 

ComparSe 4 l'imagerie spectro-temporelle de 
20 fluorescence, ce bref calcul illustratif d^montre que 
l'imagerie spectro-temporelle de transmission pr^sente 
un comptage sicurise plus important. 

On pr§cise que deux approches sont 
possibles dans la presente invention, a savoir une 
25 approche empirique et une generalisation de la 
spectrophotometrie de Beer-Lambert. 

Dans la premidre approche, on ne vise que 
la detection et dans la deuxiSme approche on vise la 
mesure . 

30 On peut relier les images 

spectrd-temporelles dSrivables et leurs dSriv§es 4 un 
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crit^re stmiologique, non optique, de nature 
statistique, par exemple le murissement d'un fruit, la 
qualite d'un proctdt UHT pour le lait, les proprietts 
mecaniques d'un plastique diffusant, la qualite d'une 
5 dispersion de ctramique. 

A titre purement indicatif et nullement 
limitatif, on a mis en oeuvre l 1 invention pour divers 
objets tels que des pommes de terre, des oeufs et des 
pommes. Dans chaque cas, on a obtenu une pile 
10 d' images derivables: 

- une premiere image constitute par un 
ensemble de points dans un repSre Nombre de SPE (dans 
un intervalle de longueurs d'onde, par exemple 700nm- 
725nm) en fonction de Temps (ps) et Temps macroscopique 

15 (s), 

- une deuxieme image constitute par un 
ensemble de points dans un repdre Nombre de SPE en 
fonction de Longueur d'onde (en nm) , et Temps 
macroscopique (s) , 

20 - une troisiSme image dite fondamentale, 

constitute par un ensemble de points (nombre de SPE) 
dams un repere Longueur d'onde (en nm) - Temps (en ps) . 

Considtrons maintenant la generalisation de 
la spectromttrie de Beer-Lambert. 

25 En gtntral, il est ntcessaire d'ttablir un 

moddle mathematique de la structure du milieu a sonder. 
Le cas le plus simple consiste & considerer I'objet & 
mesurer comme un module en* suspension isotrope 
homogene. Dans ce cas, de nombreux modules 

30 mathematiques (outre les mtthodes du type Monte -Carlo) 
ont dtja ett proposes et reposent surtout sur 
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1 Approximation de la diffusion et, plus generalement, 
sur la theorie du transfert radiatif . 

On precise que, dans le cas de deux phases 
non miscibles, la diffusion devient non linSaire quand 
5 la fraction volumique dSpasse quelques pour mille. 

Dans le cas lin£aire de la tres faible 
diffusion, un equivalent de la Loi de Beer-Lambert peut 
etre utilise, la turbidite etant alors consid6r£e cotnme 
une pseudo-absorbance: 
10 Log (I/I 0 )=-L/L* oil L* est la longueur de 

diffusion. 

Quand ce seuil est dSpasse, le regime de 
diffusion multiple commence. 

Pour une dilution de 2% ou inf^rieure a 2%, 

15 la diffusion simple s' applique. 

Les methodes de granulom^tries utilisent 
souvent un proc§d§ de dilution pour revenir au cas 
lineaire de la diffusion simple oH la theorie de Mie 
reste applicable. 

20 - Dans la present e invention, on propose de 

traiter des suspensions trds chargSes et color£es, qui 
autorisent des mesures in situ, sans contact et sans 
dilution. Dans ce cas, la generalisation de la 
spectrophotom^trie de Beer-Lambert s' applique vis-S-vis 

25 des absorbeurs avec les techniques d'imagerie 
proposees . 

Avec la technique de spectrophotom£trie 
impulsionnelle, on a par exemple Studie le cas de 
diverses concentrations d'un colorant dans du lait 
30 demi-£cremS UHT. La non utilisation des d€riv§es 
* partielles a soulignS la difficult^ d'une mesure 
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f iable. 

On a utilise un dispositif identique a 
celui qui est utilise en spectrophotomStrie classique 
(cuvette de 10mm, detection et collection centrSes sur 
5 une face) . On a utilise une procedure sequentielle mais 
un mode de double faisceau serait facile a mettre en 
oeuvre . 

L' image spectro-temporelle avec ses 
derivies partielles du continuum dans l'air puis celle 

10 du lait UHT demi-Screme non dilue sans ajout de 
colorant ont aussi 6t& enregistr^es . Puis on a ajout£ 
de faibles quantites de colorant pr^alablement dilutes 
dans l'eau. La detection de 50nanomoles/l de ce 
colorant a- et6 obtenue. 

15 Notons qu ! en spectrophotom^trie classique 

le spectre d'un colorant est tou jours enregistre par 
rapport 3. un corps de reference qui est g£n6ralement le 
solvant pur. Le lait peut etre considers ici comrae le 
solvant de reference. Le domaine spectral est choisi 

20 par rapport au pic de ce colorant mais peut §tre 
facilement transposable dans 1 'ultraviolet ou le proche 
inf rarouge . 

On a trouvS que le lait demi-6creme UHT 
diffuse beaucoup et absorbe dans le domaine bleu- vert. 
25 La comparaison avec l ! eau montre une absorption 
importante dans le domaine bleu-vert, qui est en partie 
liee aux flavines libres. 

L' invention permet de mesurer ainsi la 
moyenne de la concentration de ces colorants endogdnes 
30 sur un volume de quelques cm 3 dans l'exemple proposS. 
Mais des volumes plus importants (quelques dm 3 ) peuvent 
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§tre sondes avec une source assez puissante ou un temps 
de mesure assez long. 

L'accSs a une moyenne volumique de ce type 
de contr61e in situ est susceptible de constituer un 
5 bon index de quality du lait. 

La spectrophotom£trie impulsionnelle et 
dif f €rentielle utilis§e dans l 1 invention permet de 
mesurer une concentration d'absorbeur dans un milieu 
non limpide. Elle permet aussi de mesurer des 
10 concentrations de diffuseurs. En outre, elle permet de 
quantifier la quality du sens de la moyenne de ces 
concentrations (formation d'une moyenne ou non 
homogeneisation) . 

II convient en outre de noter que 
15 1' invention permet de realiser une imagerie spectro- 
temporelle derivable avec une seule impulsion laser- 
blanc, ce qui est impossible avec un systeme accordable 
a balayage. 

En outre, 1* invention presente l'avantage 
20 d'etre vine technique sans contact, permettant Tine 
mesure globale in situ et en ligne. Par exemple, la 
mesure en ligne de pots de yaourt ou d'objets 
diffusants peut etre realis^e dans leur propre 
contenant . 

25 L 1 invention permet en outre d'acqu£rir des 

informations sur le volume de chaque objet (et sur des 
sous-volumes absorbants, par exemple des morceaux de 
fruits), sur la qualite de la fermentation et sur la 
quality visuelle. 

30 De plus, 1' invention constitue une 

technique rapide et fiable dans le cas d'une approche 
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empirique. 

En fait, elle s'adresse de fagon generale a 
la mesure des objets diffusants. 

Dans le cas de mesures precises de 
5 concentrations d'absorbeurs (et/ou de diffuseurs) , 
1' invention peut etre consideree comme une 
generalisation de la spectrophotometrie de Beer-Lambert 
des milieux non-limpides, mais sans balayage et en mode 
impulsionnel, la spectrophotometrie de 1' invention 

10 etant en outre dif ferentielle. 

De plus, pour la mise en oeuvre de 
l 1 invention, il convient d 1 adapter la longueur impact- 
detection en fonction de 1' importance de la diffusion 
et de 1 1 absorption. 

15 Cependant, cette adaptation n f est pas 

necessaire si l'on realise des d£placements de la 
source de lumidre ou du collecteur. 

On montre dans ce qui suit les intirets des 
d£rivees partielles de 1' image spectro-temporelle que 

20 l'on utilise dans 1' invention. 

Les publications sur le transport radiatif 
en optique sont trds nombreuses [Ishimaru, 1978; 
Ishimaru A. , "Wave propagation and scattering in random 
media » NewYork, 1978.], [Mottin,2001 ; Mottin S. et 

25 Laporte P., Systdmes femtosecondes pour 1' etude des 
milieux h£t£rog£nes compartimentSs absorbants et 
diffusants dans w syst§mes femtosecondes"- Ed. P. 
Laporte, F. Salin et S. Mottin, Puse, pp. 295-310, 
2001]. Suivant 1' approximation de la diffusion bien 

30 connue, les interactions peuvent se r6duire & 
1' Equation suivante: 
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[|-(0V 2 +<^>j(ii(r,O)=5o(r,O 

5 u(r,t) est la density de photons (nombre de photons/ram 3 
ou en J/mm 3 ) . 

Le milieu sond6 se rSduit £ 2 grandeurs 

- <1/Z a >=<c[la>- Cette grandeur est homog£ne 
k une cinetique du premier ordre (en ps" 1 ) . 2a est 

10 homogdne a un temps, c est la constante universelle de 
la vitesse de la lumidre dans le vide (0,3mm/ps), 
divis£e par l'indice moyen du milieu (n=l,4). Et |^ est 
le coefficient d' absorption en mm- 1 . 

- <D> est le coefficient de dif fusivit£ 
15 optique en mm 2 /ps avec D=cD. D= (3 (lia+^g' ) ) " 1 . |V est le 

coefficient de diffusion reduit en mm" 1 . 

Cette approximation de la diffusion n'est 
utilisable que si l'on est loin des temps courts 
(c'est-a-dire si t est tr£s superieur & (cD' 1 )" 1 et si 

20 <\i a > est trds infSrieur & <|V>. 

Pour un cas simple (milieu semi-infini, <(Ha> 
et <Hs'> homogeneisables constants, un point source 
consider^ comme une distribution de Dirac spatiale, un 
point de collection considere comme une distribution de 

25 Dirac spatiale, et p la distance entre ces deux 
points), l'6quipe de Patterson [Kienle, 1997 ; Kienle 
A., Patterson M. S., Improved solutions of the steady 
state and the time-resolved diffusion equations for 
reflectance from a semi-infinite turbid medium, J. Opt. 

30 Soc. Am. A 14 246-254, 1997] a public la solution 
analytique. t=0 est 1' instant de 1' injection dans le 
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volume a partir du point source) : 



5 

avec ri=z 0/ z 0 =(jls / )" 1 /r-i=z 0 +2z b /Zb =2D (1-R d ) et Ra = 

l,440n~ 2 +0,710 n- 1 +0,668 + 0,636 n, n gtant l'indice de 
refraction qui depend de la longueur d'onde. U(t) est 
la fonction d'Heaviside. Le temps de Fourier est note 
10 t P (p,D) - p 2 /(4<D>) . 

Pour un p donnS, R(p,t) est une fonction du 
temps de vol et de la longueur d'onde. 

On peut aussi l'approximer en ' n^gligeant le 
r61e des conditions limites par: 



15 



rr[t- f 3L]:opz[X] ^ 



5/2 



avec P constant e par rapport £ t et X, 
T a (A,)=2 a/ T F (X) , toutes les deux fonctions de la 
20 longueur d'onde X {qui sont les deux fonctions 
recherchSes) . 

Cette solution analytique est une 
approximation de la distribution de probability de 
1' image spectro-temporelle. Avec le mode comptage, un 
25 ajout d'ajustement de nature poissonienne ll cette 
solution analytique est possible. 

Un des axes des caract£res innovants de la 
pr£sente invention est 1' utilisation des d£riv€es 
partielles en longueur d'onde et en temps de cette 
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image. 

Pour les cas de spectrometrie statique 
1' utilisation de la derivSe (dans ces cas, non 
partielle) est bien connue notamment pour la 
5 chimiomgtrie dans les milieux limpides. Dans la 
presente invention, dans le cadre d'une spectrometrie 
impulsionnelle (femto, pico ou nanoseconde) 
l'utilisation conjointe des derivees partielles en 
longueur d'onde et en temps de vol est nScessaire pour 
10 les milieux non limpides. 

A) L'int6rSt des derivations nidme en 
fonction du temps de vol apparait clairement ci-apres 
(on ne detaille que les derivees premiere et seconde) : 
1) Le ratio de la d£rivSe partielle en 
15 fonction du temps de vol par la fonction elle-meme est 
£gal a: 

§j(™ a & _ 1 5 | TF\X] 
image " Ta\x\ It T 2 

20 2) Le ratio de la dirivee partielle seconde 

en fonction du temps de vol par la fonction est egal k: 

1 ] 5t-2TF[Z] v 35t 2 -2ZtTF[Jl]+4TF[X\ 2 
Td[A] 2 Ta[X\t 2 4t 4 

et ainsi de suite pour les d^rivSes nidmes. 
Bri&vement on peut remarquer que 

a) la fonction en t- 5/2 e - (t/TA+TP/t > a disparu 

b) la ligne suivant X des maxima r£duits 
avec la solution d'une Equation simple du 



25 



coincide 
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second degre dont la seule solution r£aliste est : 

4 ]y 25 7fl{A] J 

5 

c) les lignes suivant X des points 
d' inflexion rgduits correspondent aux solutions d'une 
equation du quatridme degrS: 

Cette Equation a toujours 4 solutions dont 
10 deux sont reelles et positives. Ces deux derniere 
solutions rSalistes correspondent aux 2 points 
d' inflexion pour une seule longueur d'onde donnee. 

B) Pour les dSriv^es partielles eii fonction 
de la longueur d'onde, le probleme est plus dSlicat car 
15 on ne connalt pas a priori les fonctions Ta [X], TF [X] et 
Z[X]. 

1) Neanmoins Z [X] et TF[A,] sont souvent 
monotones car |V [X] 1'est tres souvent (avec 
Z[X]=2//V [X] et TF[X]~3/4p 2 fl g ' [X]/c[A]) done leurs 

20 d^rivees partielles en fonction de X sont de meme signe 
sur la fenetre spectrale. 

2) Par contre l/Ta[X] (Sgal a \i*[X]c[X]) 
presente tous les cas'de figure des fonctions continues 
avec maxima et minima. Pour la fonction \la[X) qui 

25 correspond a la density optique lin^aire de la 
spectrora§trie d' absorption des milieux limpides, vine 
litterature abondante a deja traitg 1'interSt des 
derivations niSmes (derivations non partielles) et les 
mSthodes d' analyse d6di§es. Ces approches sont limitSes 
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a un vecteur spectral a une dimension et jamais 
etendues a une image spectro-temporelle i deux 
dimensions avec un temps de vol comme seconde 
coordonn€e. L'une des activitSs innovantes ici 
5 clairement et analytiquement pr^cisSe, Slargit ces 
methodes aux derivees partielles dans le cas d'une 
imagerie d' absorption spectro-temporelle des milieux 
non limpides. 

C) De plus, un exemple d' utilisation d'une 
10 derivSe partielle croisSe est la derivation suivant X 
consecutive I la derivation suivant le temps de vol. 
Cette d£riv£e s'annule pour les lignes des maxima 
reduits • 



15 



d_ 

dX 



—(image) 



image 



= 0 



On obtient une ecriture trds simple; 



dX 



4 dX{ c[X] ) 



20 



Si l'on neglige les variations de la 
vitesse moyenne de propagation en fonction de X alors 
c[A,]=c. 

Ainsi pour les lignes des. maxima reduits 
25 1' importance de la grandeur adimensionnalisee 
Q=3/4p 2 /(ct) 2 est claire : 
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Dans toutes ces m^thodes, 1' image spectro- 
temporelle n'est considSrSe que comme un moment d'ordre 
5 z€ro II oil tous les auteurs ne la manipulent (par 
exemple passage au Log (image)) qu'en elle-m§me sans 
opgrateurs de derivation partielle. Ce sont les 
d§rivges partielles qui revelent le mieux la nature de 
cette imagerie. 

10 Ainsi nous venons de voir un exemple de 

mSthode utilisant les rapports entre les derivges 
partielles et 1' image elle-m§me. Finalement 
1' extraction de la donnee recherchee, la fonction 
llaftj* ou quelques points except ionnels (maxima, 

15 minima, zone de pente maximum, points d' inflexion. ) 
d'un milieu non limpide' s'opdre via les proc€d£s 
proposes dans la presente invention. Si l'6tude du 
milieu est plus physique (recherche de Jig' [X] et points 
exceptionnels) alors les approches sont aussi du meme 

20 ordre. Si le but est l'6tude complete du milieu non 
limpide alors les proc£d6s de 1' imagerie spectro- 
temporelle dif f £rentielle permettent aussi de qualifier 
et quantifier en m§me temps et \x B ' [X] et les 

points exceptionnels. 

25 On vient done de montrer, dans un langage 

math^matique approprie, les interSts des d§rivees 
partielles de 1' image spectro-temporelle. 

II convient de noter que pour la mise en 
oeuvre de la presente invention, on peut utiliser une ou 
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plusieurs implusions lumineuses monochromatiques, pour 
gclairer le milieu non limpide, et acquerir au moins 
une image spectro-temporelle de transmission a partir 
de ce milieu, en rSalisant un balayage spectral. 
5 Mais on peut aussi utiliser une ou 

plusieurs impulsions lumineuses S large spectre, pour 
eclairer le milieu non limpide, et acquerir au moins 
une image spectro-temporelle de transmission a partir 
de ce milieu, sans balayage spectral. 

10 Par « impulsion lumineuse k large 

spectre x>, on entend une impulsion lumineuse dont le 
spectre est superieur ou 6gal k la fenStre spectrale 
des moyens d' acquisition, a savoir le polychromateur et 
la camera & balayage de fente dans un exemple donne 

15 plus haut. 
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REVENDICATIONS 

1. Proc£d6 de spectrophotomStrie 
diffSrentielle, pour analyser un milieu non limpide 
(11) , ce proc£d6 etant caract6ris6 en ce que : 

5 - on eclaire le milieu non limpide par au 

moins une impulsion lumineuse permettant 1' utilisation 
ultSrieure de la d^rivSe partielle, en fonction de la 
longueur d'onde, d'au moins une image spectro- 
temporelle de transmission acquise a partir du milieu 

10 ainsi Sclaire, 

- on acquiert, a travers au moins un 
collecteur de lumiere, a partir du milieu ainsi 
Eclaire, au moins une image spectro- temporelle de 
transmission, en mode comptage permettant 1' utilisation 

15 ultSrieure des derivies partielles de 1' image en 
fonction de la longueur d'onde et du temps de vol de 
1' impulsion lumineuse, et 

on traite 1' image et ses derivees 
partielles en fonction de la longueur d'onde et du 

20 . temps de vol pour acqugrir des informations sur le 
milieu non limpide. 

2. Procede selon la revendication 1, dans 
lequel on exploite, en vue de qualifier le degr6 
d'homog£n6isation du milieu sonde (11) ou de detecter 

25 une singularity d' absorption et/ou de diffusion, les 
derivees partielles liees aux variations en temps de 
vol et en spectre des photons diffusifs et ergodiques 
non absorbes sur un large domaine spectral par le 
milieu, ces photons diffusifs et ergodiques Stant 6mis 

30 par le milieu lors de 1 »€clairement de celui-ci. 
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3. Proc€d£ selon la revendication 1, dans 
lequel on utilise conjointement les dSrivees partielles 
liees au temps de vol de l r impulsion lumineuse et aux 
donnSes spectrales pour Stablir une carte d'identite 

5 spectro-temporelle du milieu non limpide (11) . 

4. ProcSdS selon la revendication 1, dans 
lequel on §claire le collecteur sans sonder le milieu 
(11) et en meme temps le milieu non limpide par une ou 
plusieurs impulsions lumineuses autorisant une imagerie 

10 spectro-temporelle derivable par rapport & la longueur 
d'onde et au temps, deux pics simulant un double 
faisceau. 

5. Dispositif de spectrophotometrie 
dif ferentielle, pour analyser un milieu non limpide 

15 (11), ce dispositif etant caract#ris6 en ce qu f il 
comprend : 

- une source lumineuse impulsionnelle (8) 
pour eclairer le milieu non limpide et permettant 
1' utilisation ult§rieure de la derivee partielle, en 

20 fonction de la longueur d'onde, d'une image spectrale 
et temporelle de transmission acguise k partir du 
milieu ainsi 6clair£, 

- des moyens (18) d 1 acquisition, a partir 
du milieu ainsi SclairS, d f une image spectrale et 

25 temporelle de transmission, en mode de comptage, 
permettant 1' utilisation ult6rieure des d§rivies 
partielles de 1' image en fonction de la longueur d'onde 
et du temps de vol del' impulsion lumineuse, et 

- des moyens (20) de traitement de cette 
30 image, considSrSe comme le moment d'ordre zero, et de 

ses derivees partielles en fonction de la longueur 
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d'onde et du temps, pour acquerir des informations sur 
ce milieu. 

6. Dispositif selon la revendication 5, 
dans lequel la source lumineuse impulsionnelle (8) 

5 comprend des moyens de g€n6ration non lineaire 
d 1 impulsions lumineuses femtosecondes ou picosecondes 
permettant 1' utilisation de la dSrivSe partielle en 
fonction de la longueur d'onde. 

7. Dispositif selon la revendication 5, 
10 dans lequel la source lumineuse impulsionnelle comprend 

des moyens (8) de ggnSration et d' amplification de 
continuum femtoseconde ou picoseconde ou nanoseconde 
d f un continuum, le continuum permettant directement 
l'utilisation de la deriv£e partielle en fonction de la 
15 longueur d'onde. 

8. Dispositif selon la revendication 5, 
dans lequel les moyens d 1 acquisition (18) comprennent 
une camera k balayage de fente en mode comptage. 

9. Dispositif selon la revendication 5, 
20 dans lequel les moyens d r acquisition comprennent une 

camera (18) & balayage de fente de comptage de photo- 
electron unique en mode d' operation monocoup ou en mode 
d' operation synchro-scan ou en mode analogique. 

10. Proc6d£ de spectrophotomStrie 
25 dif f^rentielle selon la revendication 1, permettant 

l'etablissement d'une veritable carte d' identity 
optique du volume sonde, c'est-S-dire d'une fiddle 
signature de la nature statistique du milieu (11) vis- 
£-vis des teneurs plus ou moins homogSnes des 
30 diffuseurs et des absorbeurs, cette carte d'identite 
optique 6tant sous la forme d'une ou plusieurs images 
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spectro-temporelles permettant d' avoir accds 
simultangment aux distributions teraporelles pour une 
fen§tre spectrale donn§e, aux distributions spectrales 
pour un intervalle de temps donne, aux derivSes 
5 partielles de ces deux distributions ainsi qu'S. leurs 
int€grales. 

11. ProcSde de spectrophotom£trie et 
tomographie selon la revendication 1, dans lequel on 
utilise conjointement les d6rivees partielles, par 

10 rapport k la longueur d'onde, au temps de vol et k 
l'espace, des images spectro-temporelles et des modes 
monopoints avec balayage spatial ou des modes 
multipoints commutes des zones d' injection et/ou des 
zones de collection de la lumiere, l'utilisation des 

15 derivees partielles de type espace 6tant alors possible 
et permettant de traiter certains cas de density non 
constante dans le volume sonde, ce proc§d§ de. 
tomographie avec des dif f 6rentielles d' image spectro- 
temporelle permettant, d'une part, un rep§rage d'une 

20 singularity de la concentration d'absorbeurs et/ou de 
diffuseurs et, d' autre part, une identification 
raoleculaire de ces absorbeurs, les injections et 
collections de la lumidre pouvant §tre realisees soit & 
la surface soit au sein du volume, 

25 12. Proced£ de spectrophotomStrie 

diffgrentielle selon la revendication 1, dans lequel on 
utilise conjointement les derivees partielles, par 
rapport k la longueur d'onde et au temps, de l'imagerie 
spectro-temporelle de transmittance et le mode de 

30 comptage par binarisation de 1' image puis detection 
d'une zone de pixel attribuable k un photo- Electron et 
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reduction de cette zone k un seul pixel allume, ou a 
une echelle subpixel, afin d'augmenter la dynamique et 
de qualifier la mesure en photo-electron unique, 

13. Proc6de de spectrophotom^trie selon la 
5 revendication 1, dans lequel on exploite, en vue de 

mesurer de trds faibles variations d'un absorbeur, les 
derivges partielles en temps et en longueur d'onde de 
la queue temporelle li£e aux photons ergodiques 
diffuses par le milieu (11) lors de 1 1 eclairement de 
10 celui-ci. 

14. Procede de spectrophotom€trie 
different ielle selon la revendication 1, dans lequel on 
exploite, en vue de mesurer ou d'homog£n§iser des 
teneurs d' absorbeurs et/ou dif f useurs , 1 ' imagerie 

15 spectro- temporelle et les operateurs d/dt, 9*/dt 2 , d/dX, 
d*/dX 2 , cf*/dtdX jusqu'aux ordres superieurs, ces 
operateurs s'appliquant k 1' image spectro- temporelle. 

15 . ProcSde de spectrophotom£trie 
dif fSrentielle selon la revendication 1, dans lequel on 

20 utilise conjointement les fluctuations des images 
spectro-temporelles de transmittance et les operateurs 
associis afin d'effectuer des mesures dynamiques 
d'opacimetrie, de colorimetrie et de granulomgtrie . 



25 
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